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PSEUDO-NOISE (P-N) SEQUENCE TEST 


DECUS Program Library Write-up 


DEC US No. 8-304 


ABSTRACT 

This program can be used to determine the statistical characteristics of a pseudo-random se¬ 
quence generator. The probability density function can be studied with the aid of an ampli¬ 
tude histogram, and measures of the "randomness" can be determined with correlation dia¬ 
grams, sample function displays, and scattergrams. All diagrams are displayed on the CRT. 

REQUIREMENTS 

Storage 


Program: 0-1,16-171,200-356,400-572,600-713 
Date: 2000-5777* 

Axis Points: 1000-1777* 

*Data and Axis Points are stored as (x,y) co-ordinates, ( x 0 yo x l y l.) 

Equipment 

PDP-8 with EAE, display, A/D converter, and external interrupt. The 6424 instruction seen 
in this program resets the clock interrupt at the Carleton installation. It is used to time 
sampling of the pseudo-random sequences and should be replaced by the appropriate instruction 
or NOP when implemented at a different installation. 

USAGE 

Calling Sequence 

The program is called by setting the switch register to 0200 and pressing "load address. " The 
switches are reset for the desired test configuration and "start" is pressed. 

Switch Settings 

(I) The type of test to be used is determined by switches 1 and 2. 


0 Amplitude Histogram 

1 Correlation Diagram 

0 Scattergram 

1 Sample Function Display 











(II) The input mode is determined by switch 0. 


0 Software* 

1 A/D converter (channel 0) 

*Unless changes are made, the program will extract pseudo-random numbers from its own 
reference software generator, RAND, stored in locations 0136-0156. It is a mixed congru- 
ential generator whose sequence is defined by the following algorithm. 

x. + j = x.(2^+3) + 3 modulo 2^ 

where Xq*1 . 

If a different software generator is required, it must be read into unused locations in memory 
and a JMS to its starting address inserted in location 0074. 

(III) When the correlation diagram is used, the shift (minus one) between the samples plotted 
as co-ordinate pairs is determined by switches 3-11. Hence, the maximum shift is 512. 

(IV) A change in the test being run can be effected by resetting the switch register for the 
new test, and depressing any teletype key. 

RESTRICTIONS 

Allowable Range 

This program is intended for testing 10 bit pseudo-random numbers, considered to be in the 
low order 10 bits of the accumulator. 

DESCRIPTION 

Discussion 

(I) Amplitude Histogram - This is a free running accumulation of the number of occurrences 
of a particular number in a finite record of a pseudo-random sequence. The maximum accumu¬ 
lation is determined by the rise increment, IN. Each time a particular number occurs, its 
accumulation counter is incremented by IN. With IN set at 10 “8.-., the maximum accumu- 

o IU 

lotion is 177g = 128jq. A visual indication of the first order probability density function of the 
sequence under is given by this display. 

For pseudo-random number generators with uniform distributions, the mean accumulation and 
variance (on the mean accumulation) may be determined as follows: 

Stop the histogram test when a suitable display has been obtained (an approximate mean 
accumulation more than half way up the vertical axis). The mean accumulation will appear 
in the ACC, as defined by: 





1024 


y =1/1024 
'n 


i=l 


where y. is the value of the i** 1 accumulation counter. This must be normalized by N, the 

number of samples, to be compared with the theoretical estimate (assuming N independent 
trials of the generator). 


y f ~ 1/1024 (for 1024 quantization levels) 


Press "continue." The variance on the mean accumulation will appear, high order 12 bits 
in the ACC, low order 12 bits in the MQ: 


1024 


2 

c r = 


i=1 


(y,-yf 


2 

This must be divided by (1024*N ) for proper normalization. It may then be compared with 
the theoretical estimate of the variance on the mean accumulation (again assuming N in¬ 
dependent trials). 


<T = 1/(1024* N) 


Press continue. " The number of samples taken for the histogram will appear, high order 
12 bits in the ACC, low order 12 bits in the MQ. 


(II) Correlation Diagram - A visual indication of the correlation between successive samples 
may be obtained from the correlation diagram. It is a display of (1024+(s+l)) samples. Co¬ 
ordinates are formed from (r.,r. + ^j) where the r. are samples of the pseudo-random sequence. 


(Ill) Scattergram — A measure of the "randomness" of the sequence under test may be obtained 
from the scattergram or successive pairs diagram. Co-ordinates are formed by plotting suc¬ 
cessive pairs of samples of the sequence as (x,y) co-ordinates as follows: 


( r 0 7 r l)' ( r 2 ,l V / ( r 4 ,r 5^ 


This display is free-running and takes in 1024 such pairs before overwriting the oldest pairs. 
Further details of the use of this diagram are given by Chambers (1). 

(IV) Sample Function Display - Further indications of the structure or patterns in the sequence 
are provided by the display of a finite record of the output of the generator under test. In 
this display, (i,r.) is plotted as an (x,y) co-ordinate, where r. is a sample of the pseudo¬ 
random sequence, and i is the sample index. The display is 1 free-running and takes in 1024 
samples before overwriting the oldest samples. 
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EXECUTION TIME 


The histogram test has an approximate maximum frequency of 6 khz. The other tests have 
approximate maximum interrupt frequencies of 4 khz. 

These frequencies are quoted for sampling from an A/D converter with a 50 usee, conversion 
time. 

REFERENCES 

1. R. P. Chambers, "Random-number Generation On Digital Computers, " IEEE Spectrum, 
February 1967, pp. 48-56. 
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/PSEUDO-RANDOM SEQUENCE GENERATOR 
/AMPLITUDE HISTOGRAM, CORRELATION 
/DIAGRAM (SCATTERGRAM), AND SAMPLE 
/FUNCTION OBSERVATION. 

*0000 


0000 

0000 

0000 

0001 

5100 

JMP SER / JUMP TO INTERRUPT SERVICE 
*0016 

0016 

0000 

MEM,0 

0017 

0000 

MEMI,0 

0020 

0000 

N0,0 

0021 

0000 

N01,0 

0022 

0000 

XMIN,0 

0023 

0000 

MTEMP,0 

- 0024 

0000 

X,0 

V 0025 

0000 

Y,0 

0026 

0400 

THOU,040,0/ NUMBER OF AXIS POINTS 

0027 

7600 

CNT, 7600/ NO. OF POINTS ON 1 AXIS 

0030 

0000 

CNT1,0 

0031 

0010 

INC,0010/ INCREMENT BETWEEN POINTS 

0032 

0000 

A0,0/ ORIGIN 

0033 

0000 

B0,0/SETTINGS 

0034 

0000 

AC,0/ ACC-TEMPORARY STORAGE 

0035 

0000 

PL1,0000 

0036 

0200 

PL2,0S00 

0037 

0777 

M0/0777 / INITIAL AXIS ADDRESS MINUS ONE 

0040 

2000 

X0,2000/ INITIAL X ADDRESS 

0041 

2001 

Y0,2001 / INITIAL Y ADDRESS 

0042 

0000 

XLOC,0 

0043 

0000 

YLOC,0 

0044 

0000 

CNTR,0 

0045 

2000 

NRAND,2000 

W 0046 

0000 

IR,0/ SAMPLE INDEX 

0047 

0000 

S,0/ SEPARATION MINUS ONE 

0050 

0000 

T1,0 

0051 

1777 

MASK, 1777 

0052 

0777 

MASKS ,0777 

0053 

0177 

MAS KV, 0177 

0054 

0003 

MAS KH,0003 

0055 

0535 

EXEC, AXIS 

0056 

1000 

CORR, 1000/ Y AXIS SHIFT ON HISTOGRAM 

0057 

0336 

HT, HSTINT 

0060 

0400 

SCI,SCAT1 / INTERRUPT SERVICE POINTERS 

0061 

0417 

SC2, SCAT2 

0062 

0442 

SC3, SCAT3 

0063 

0456 

SC4* WAIT 

0064 

0460 

SCP, SEQ 

0065 

0507 

SA1,SAMP1 

0066 

7777 

SOFT, 7777 

0067 

0000 

HIGH,0000 
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0070 0000 
0071 0000 

0072 0000 
0073 4157 

0074 4136 

0075 0000 

0076 0000 
0077 0000 
0100 3034 

0101 6424 

0102 6031 

0103 74X0 

0104 5436 

0105 4435 

0106 7300 

0107 1034 

0110 6001 
011 1 5400 

0112 6077 
0113 6424 

0114 6001 

0115 7300 

0116 1016 
0117 3023 

0120 1020 
0121 7041 

0122 3021 

0123 1416 

0124 6052 

0125 7300 

0126 1416 

0127 6066 

013^ 7300 

0131 2021 

0132 5123 

0133 1023 

0134 3016 

0135 5115 

0136 0000 

0137 7621 

0140 1155 

0141 7413 

0142 0005 
0143 1155 

0144 1155 

0145 1155 

0146 1156 

0147 3155 

0150 1155 

0151 7012 

0152 0051 
0153 5536 


LOW,0000 
H OR,0000 
LOR,0000 

CHGE, JMS A2D / SOFTWARE,/\/D CONVERTER 
SET, JMS RAND / SAMPLING LINKS 
CNIX,0 

NPT1 ,0 / SAMPLES COUNTER, HIGH ORDER 
NPT2,J0 / LOW ORDER 

SER, DCA AC / INTERRUPT SERVICE LINK 
6424 / RESET CLOCK 
6031 / KEYBOARD INTERRUPT? 

SKP 

JMP I PL2 / YES, SET UP NEW TEST 
JMS I PL1 / NO, SERVICE INTERRUPT 
CLACLL 

TAD AC / RESTORE ACC 
ION / ENABLE INTERRUPT 
JMP I 0000 / CONTINUE DISPLAY 
TRACE, 6J077 / DISPLAY ROUTINE-SET INTENSITY 
6424 / RESET CLOCK 
ION / ENABLE INTERRUPT 
START, CLA CLL 

TAD MEM / PUT INITIAL ADDRESS MINUS 
DCA MTEMP / ONE IN TEMP STORE 
TAD NJ0 

CIA / SET POINTS COUNTER 
DCA Njefl 

C, TAD I MEM 

6052 
CLA CLL 

TAD I MEM / DISPLAY CO-ORDINATE 
6066 
CLA CLL 

ISZ N01 / FINISHED? 

JMP C / NO, CONTINUE 
TAD MTEMP / YES, REINITIALIZE 
DCA MEM 

JMP START / RESTART DISPLAY 
RAND, 0 / MIXED CONGRUENTIAL 

CLA MQL / PSEUDO-RANDOM NUMBER 

TAD RN / GENERATOR 

SHL / R(l+l)=R(l)*(2T6+3)+3 

0005 / MODULO 2T12 

TAD RN 

TAD RN 

TAD RN 

TAD THREE 

DCA RN 

TAD RN 

RTR 

AND MASK / TAKE HIGH ORDER 10 BITS 
JMP I RAND / FOR DISPLAY TESTS 




6 



0154 

0001 

0155 

0000 

0156 

0003 

0157 

0000 

0160 

• 7300 

0161 

6542 

0162 

6532 

0163 

6531 

0164 

5163 

0165 

6534- 

0166 

5557 

0167 

0000 

0170 

4136 

0171 

5567 

# 0200 

6032 

0201 

7300 

0202 

1154 

0203 

3155 

0204 

1074 

0205 

3170 

0206 

1037 

0207 

3016 

0210 

1045 

0211 

1026 

0212 

3020 

0213 

7604 

0214 

7104 

0215 

7620 

0216 

5221 

0217 

1073 

g- 0220 

3170 

v 0221 

7604 

0222 

7006 

0223 

7006 

0224 

0054 

0225 

7450 

0226 

5236 

0227 

1066 

0230 

7450 

0231 

5244 

0232 

1066 

0233 

7650 

0234 

5257 

0235 

5265 

0236 

1056 

0237 

3022 

0240 

3024 

0241 

1057 


RN|0,0001 
RN,0000 
THREE, 0002 

A2D, 0 / SAMPLE SEQUENCE VIA A/D CONVERTER 
CLA CLL 

ADSC / SET TO CHANNEL 0 

ADCV 

ADSF 

JMP .-1 / SAMPLE SEQUENCE 
ADRB 

JMP I A2D / RETURN 
FETCH, 0 / FETCH A SAMPLE 
Gl, JMS RAND 

JMP I FETCH / RETURN 
*0200 / START OF TESTS 
6032 / CLEAR KEYBOARD FLAG 
CLA CLL 

TADRN0 / INITIALIZE SOFTWARE 
DCA RN / GENERATOR 
TAD SET 

DCA Gl / SET SOFTWARE-A/D CONVERTER CHOICE 
TAD MQ / SET INITIAL DATA STORAGE ADDRESS 
DCA MEM 

TAD NRAND / SET POINTS COUNTER 
TAD THOU 
DCA N0 
LAS 

CLL RAL 

SNL CLA / READ BIT0; CHOOSE INPUT MODE 
JMP TSTRD / SOFTWARE REQUIRED (0) 

TAD CHGE / A/D CONVERTER REQUIRED (1) 

DCA Gl 
TSTRD, LAS 
RTL 
RTL 

AND MASKH / READ BITS 1,2; CHOOSE TEST 
SNA 

JMP HIS / HISTOGRAM REQUIRED (00) 

TAD SOFT 
SNA 

JMP SCA1 / SCATTERGRAM REQUIRED (01) 

TAD SOFT 
SNA CLA 

JMP SCAP / SUCCESSIVE PAIRS PLOT (02) 

JMP SCA2 / SAMPLE FUNCTION REQUIRED (03) 

HIS, TAD CORR / HISTOGRAM PRELIMINARY 
DCA XMIN / CORRECT Y AXIS MARKER 
DCA X / CLEAR AXIS COUNTER 
TAD HT / SET UP HISTOGRAM INTERRUPT 


7 



,0242 

3035 


DCA PL1 / SERVICE LINK 

0243 

5316 


JMP HISTGM / PROCEED WITH HISTOGRAM 

0244 

4272 

SCA1, 

JMS CLEAR / SCATTERGRAM PRELIMINARY 

0245 

1060 


TAD SCI / SET UP SCATTERGRAM INTERRUPT 

0246 

3035 


DCA PL1 / SERVICE LINK 

0247 

7604 


LAS / DETERMINATION SEPARATION (S) 

0250 

0052 


AND MASKS / FROM BITS 3-11 

0251 

7001 


IAC 

0252 

3047 


DCA S 

0253 

1047 


TAD S 

0254 

7041 


CIA / SET HISTORY COUNTER 

0255 

3044 


DCA CNTR 

0256 

5455 


JMP 1 EXEC / JUMP TO AXIS GENERATION 

0257 

4272 

SCAP, 

JMS CLEAR / SUCCESSIVE PAIRS PLOT 

0260 

1064 


TAD SCP / PRELIMINARY 

0261 

3035 


DCA PL1 

0262 

7240 


CLA CMA / SET XY COUNTER 

0263 

3075 


DCA CNTX 

0264 

5455 


JMP 1 EXEC / JUMP TO AXIS GENERATION 

0265 

4272 

SCA2, 

JMS CLEAR / SAMPLE FUNCTION PRELIMINARY 

0266 

1065 


TAD SA1 / SET UP SAMPLE FUNCTION INTERRUPT 

0267 

3035 


DCA PL1 / SERVICE LINK 

0270 

3046 


DCA IR / CLEAR SAMPLE INDEX 

0271 

5455 


JMP 1 EXEC / JUMP TO AXIS GENERATION 

0272 

0000 

CLEAR, 

0 / INITIALIZATION 

0273 

1045 


TAD NRAND 

0274 

7104 


CLL RAL 

0275 

7041 


CIA 

0276 

3021 


DCA N01 

0277 

1040 


TAD X0 

0300 

3042 


DCA XLOC 

0301 

3442 


DCA 1 XLOC / CLEAR DATA REGISTERS 

0302 

2042 


ISZ XLOC 

0303 

2021 


ISZ N01 

0304 

5301 


JMP .-3 

0305 

1040 


TAD X0 / CLEAR X ADDRESS 

0306 

3042 


DCA XLOC 

0307 

1041 


TAD Y0 

0310 

3043 


DCA YLOC / CLEAR Y ADDRESS 

0311 

3022 


DCA XMIN / CLEAR Y AXIS MARKER 

0312 

1045 


TAD NRAND 

0313 

7041 


CIA / SET POINTS COUNTER 

0314 

3050 


DCA T1 

0315 

5672 


JMP 1 CLEAR / RETURN 

0316 

1045 

HISTGM, 

TAD NRAND / HISTOGRAM SETUP 

0317 

7041 


CIA / SET POINTS COUNTER (1024) 

0320 

3021 


DCA N01 

0321 

1040 


TAD X0 

0322 

3042 


DCA XLOC / SET REGISTERS 




8 



$323 
$324 
$325 
$326 
$327 
$33$ 
$331 
$332 
$333 
$334 
$335 
$336 
$337 
$34$ 
$341 
' $342 
$343 
$344 
$345 
$346 
$347 
$35$ 
$351 
$352 
$353 
$354 
$355 
$356 


1 $24- 
3442 
3042 
3442 
2042 
2024 
2$21 
5323 
3$76 
3$77 
5455 

$$$$ 

73$$ 

1 $77 
7$$1 
3$77 
7004 
1 $76 
3$76 
4167 
71 $4 
1$41 
3$43 
1443 
1356 
3443 
5736 
$$ 1 $ 


CLR, 


HSTINT, 


IN,$01$ 


0400 

$$$$ 

SCAT1 

0401 

2044 


4,$4$2 

7410 


^ $403 

5211 


$4$4 

4167 

FI, 

0405 

3442 


0406 

2$42 


0407 

2042 


0410 

56$$ 


0411 

1$61 

FIN, 

0412 

3$35 


0413 

1$5$ 


$414 

1$47 


$415 

3$5$ 


$416 

5204 


$417 

$$$$ 

SCAT2, 

$420 

2050 


$421 

7410 


$422 

5234 


$423 

4167 

F2, 

$424 

3443 



/ FETCH SAMPLE 

TIMES TWO FOR PROPER STORAGE 


TAD X 
DCA I XLOC 
ISZ XLOC 
DCA I XLOC 
ISZ XLOC 
ISZ X 
ISZ N$1 
JMP CLR 

DCA NPT1 / CLEAR SAMPLES COUNTER 
DCA NPT2 

JMP I EXEC / JUMP TO AXIS GENERATION 
$ / HISTOGRAM INTERRUPT SERVICE 
CLA CLL 
TAD NPT2 

I AC / INCREMENT SAMPLES COUNTER (24 BIT) 

DCA NPT2 

RAL 

TAD NPT1 
DCA NPT1 
JMS FETCH 
CLL RAL / 

TAD Y$ 

DCA YLOC 

TAD I YLOC / UPDATE DISPLAY 
TAD IN 
DCA I YLOC 

JMP I HSTINT / RETURN 
/ RISE INCREMENT OF HISTOGRAM 
*0400 

0 / SCATTERGRAM SERVICE 1 
ISZ CNTR 
SKP 

JMP FIN / TAKE IN S SAMPLES AND STORE 
JMS FETCH / IN X ARRAY 
DCA I XLOC 
ISZ XLOC 
ISZ XLOC 
JMP I SCAT1 

TAD SC2 / FINISHED. CHANGE INTERRUPT 
DCA PL1 / LINK 
TAD T1 
TAD S 
DCA T1 
JMP FI 

$ / SCATTERGRAM SERVICE 2 
ISZ T1 
SKP 

JMP FIN2 / TAKE IN (NRAND-S) 

JMP FETCH / SAMPLES AND STORE 
DCA I YLOC / IN X AND Y ARRAYS 
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$425 

1443 


TAD 1 YLOC 

$426 

3442 


DCA 1 XLOC 

$427 

2042 


ISZ XLOC 

$43$ 

2042 


ISZ XLOC 

$431 

2043 


ISZ YLOC 

$432 

2043 


ISZ YLOC 

$433 

5617 


JMP1SCAT2 

$434 

1$62 

FIN2, 

TAD SC3 / FINISHED. CHANGE INTERRUPT 

$435 

3$35 


DCA PL1 / LINK 

$436 

1$47 


TAD S 

$437 

7041 


CIA 

$44$ 

3$44 


DCA CNTR 

$441 

5223 


JMP F2 

$442 

$$$$ 

SCAT3, 

$ / SCATTERGRAM SERVICE 3 

$443 

2044 


ISZ CNTR 

$444 

741$ 


SKP 

$445 

5253 


JMP FIN3 / TAKE IN S SAMPLES AND 

$446 

4167 

F3, 

JMS FETCH / STORE IN Y ARRAY 

$447 

3443 


DCA 1 YLOC 

$45$ 

2043 


ISZ YLOC 

$451 

2043 


ISZ YLOC 

$452 

5642 


JMP 1 SCAT3 

$453 

1$63 

FIN3, 

TAD SC4 / FINISHED. CHANGE INTERRUPT 

$454 

3$35 


DCA PL1 / LINK 

$455 

5246 


JMP F3 

$456 

$$$$ 

WAIT, 

$ / INTERRUPT WAIT LOOP 

$457 

5656 


JMP 1 WAIT 

$46$ 

$$$$ 

SEQ, 

$ / SUCCESSIVE PAIRS SERVICE 

$461 

4167 


JMS FETCH / FETCH SAMPLE 

$462 

2$75 


ISZ CNTX / X OR Y 

$463 

5272 


JMP YST 

$464 

3442 


DCA 1 XLOC / STORE X 

$465 

2042 


ISZ XLOC 

$466 

2042 


ISZ XLOC 

$467 

7344 


CLA CLL CMA RAL / RESET COUNTER 

$47$ 

3075 


DCA CNTX 

$471 

566$ 


JMP 1 SEQ / RETURN 

$472 

3443 

YST, 

DCA 1 YLOC / STORE Y 

$473 

2043 


ISZ YLOC 

$474 

2043 


ISZ YLOC 

$475 

205$ 


ISZ T1 / ADDRESS RESET NEEDED? 

$476 

566$ 


JMP 1 SEQ / NO, RETURN 

$477 

104$ 


TAD X$ / YES 

$5$$ 

3042 


DCA XLOC 

$5$1 

1041 


TAD Y$ 

$5$2 

3043 


DCA YLOC 

$5$3 

1045 


TAD NRAND 

$5$4 

7041 


CIA 

$5$5 

305$ 


DCAT1 
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0506 5660 

0507 0000 

0510 7300 

051 1 2050 

0512 5323 

0513 1045 

0514 7041 

0515 3050 

0516 1040 

0517 3042 

0520 3046 

0521 1041 

0522 3043 

0523 4167 

• 0524 3443 

0525 2043 

0526 20 43 

0527 1046 

0530 3442 

0531 2042 

0532 2042 

0533 2046 

0534 5707 

0535 101& 

0536 3023 

0537 4365 

0540 1022 

0541 3416 

0542 1025 

0543 3416 

0544 1025 

(r 0545 1031 

0546 3025 

0547 2030 

0550 5340 

0551 4365 

0552 1024 

0553 3416 

0554 1024 

0555 1031 

0556 3024 

0557 3416 

0560 2030 

0561 5352 

0562 1023 

0563 3016 

0564 5112 

0565 0000 

0566 1027 


JMP I SEQ / RETURN 

SAMP1, 0 / SAMPLE FUNCTION INTERRUPT SERVICE 

CLA CLL 
ISZ T1 
JMP GG3 

TAD NRAND / RESET ADDRESSES 

CIA 

DCA T1 

TAD X0 

DCA XLOC 

DCA IR 

TAD Y0 

DCA YLOC 

GG3, JMS FETCH / FETCH SAMPLE 
DCA I YLOC 

ISZ YLOC / UPDATE DISPLAY 

ISZ YLOC 

TAD IR 

DCA I XLOC 

ISZ XLOC 

ISZ XLOC 

ISZ IR 

JMP I SAM PI / RETURN 
AXIS, TAD MEM / AXIS GENERATION 
DCA MTEMP 

JMS RSET / INITIALIZE 

A, TAD XMIN 
DCA I MEM 

TAD Y / STORE VERTICAL AXIS POINTS 
DCA I MEM 
TAD Y 

TAD INC / INCREMENT COUNTER 
DCA Y 
ISZ CNT1 
JMP A 

JMS RSET / INITIALIZE 

B, TAD X 

DCA I MEM / STORE HORIZONTAL AXIS POINTS 
TAD X 

TAD INC / INCREMENT COUNTER 

DCA X 

DCA I MEM 

ISZ CNT1 

JMP B 

TAD MTEMP 
DCA MEM 

JMP TRACE / JUMP TO DISPLAY 
RSET, 0 / INITIALIZATION ROUTINE 
TAD CNT 
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0567 

3030 


DCA CNT1 

0570 

3024 


DCA X 

057 1 

3025 


DCA Y 

0572 

5765 


JMP 1 RSET 
*0600 

0600 

4247 

MEAN, 

JMS INIT / CALCULATE SAMPLE MEAN AND VARIANCE 

0601 

1443 

GO, 

TAD 1 YLOC / ON HISTOGRAM 

0602 

7010 


RAR 

0603 

7012 


RTR / CORRECT FOR HISTOGRAM RISE INCREMENT 

0604 

0053 


AND MASKV / 0177 IS MAX ACCUMULATION (128) 

0605 

3070 


DCA LOW 

0606 

4263 


JMS SUM / SUM ACCUMULATIONS 

0607 

2043 


ISZ YLOC / AND COMPUTE MEAN 

0610 

2043 


ISZ YLOC 

0611 

2021 


ISZ N01 

0612 

5201 


JMP GO 

0613 

4304 


JMS AV / AVERAGE AND DISPLAY MEAN IN ACC 

0614 

7402 


HLT 

0615 

7041 

VAR, 

CIA / INVERT & STORE MEAN 

0616 

3024 


DCA X 

0617 

4247 


JMS INIT / INITIALIZE 

0620 

1443 

LOOP, 

TAD 1 YLOC 

0621 

7010 


RAR 

0622 

7012 


RTR / CORRECT FOR HISTOGRAM RISE INCREMENT 

0623 

0053 


AND MASKV / 0177 

0624 

1024 


TAD X 

0625 

7510 


SPA / COMPUTE (X-MEAN)T2 

0626 

7041 


CIA 

0627 

3277 


DCA SS 

0630 

1277 


TAD SS 

0631 

4275 


JMS SQ 

0632 

4263 


JMS SUM / SQUARE & SUM 

0633 

2043 


ISZ YLOC 

0634 

2043 


ISZ YLOC 

0635 

2021 


ISZ N01 

0636 

5220 


JMP LOOP 

0637 

1072 


TAD LOR 

0640 

7421 


MQL / DISPLAY 24 BIT UNAVERAGED VARIANCE 

0641 

1071 


TAD HOR / HIGH ORDER IN ACC 

0642 

7402 


HLT / LOW ORDER IN MQ 

0643 

1077 


TAD NPT2 / DISPLAY NUMBER OF 

0644 

7421 


MQL / SAMPLES (24 BITS) 

0645 

1076 


TAD NPT1 / HIGH ORDER IN ACC 

0646 

7402 


HLT / LOW ORDER IN MQ 

0647 

0000 

INIT, 

0 / INITIALIZATION ROUTINE 

0650 

7300 


CLA CLL 

0651 

3071 


DCA HOR 

0652 

3072 


DCA LOT 

0653 

3067 


DCA HIGH 

0654 

3070 


DCA LOW 
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jto i3ii 
0656 704] 

0657 3021 

0660 1041 

0661 3043 

0662 5647 

0663 0000 SUM, 

0664 7100 

0665 1070 

0666 1072 

0667 3072 

0670 7204 

0671 1071 

0672 1067 

• 0673 3071 

0674 5663 

0675 0000 SQ, 

0676 7425 

0677 0000 SS, 

0700 3067 

0701 7501 

0702 3070 

0703 5675 

0704 0000 AV, 

0705 1072 

0706 7421 

0707 1071 

0710 7407 

0711 1777 N00, 

0712 7701 

i 0713 5704 


TAD N00 
CIA 

DCA Njefl 

TAD Y0 

DCA YLOC 

JMP I I NIT 

0/24 BIT ADDER 

CLL / ENTER WITH NO. IN HIGH 

tad low 
tad lor 

DCA LOR 
GLK 

TAD HOR 
tad HIGH 
DCA HOR 
JMPI SUM 

0/24 BIT SQUARING ROUTINE 
MQLMUY 
0 


& LOW 


DCA HIGH 
MQA 
DCA LOW 

JMP I SQ / EXIT WITH PRODUCT 
0/24 BIT DIVIDER 
TAD LOR 
MQL 


IN HIGH & LOW 


TAD HOR 
DVI 


^7~77 / 1024 QUANTIZATION LEVELS 
CLA MQA 

JMP IAV / EXIT WITH 12 BIT QUOTIENT IN ACC 
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